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We investigate the role of radiative cooling within the core of large X-ray clusters 
using multi-phase hydrodynamics. We developed for that purpose a spherically 
symmetric hydrodynamical code, coupled to a "fluid model" that describes ac- 
curately the dark matter component. Cooling is included using a self-consistent 
multi-phase approach, leading to "cooled gas" mass deposition throughout the 
flow. We simulate the collapse and the subsequent evolution of a Coma-like X-ray 
cluster, avoiding the well-known "cooling catastrophe". The total mass profile of 
our simulated cluster ipivery similar to the "universal" profile proposed by Navarro, 
Frenk & White (1995)EJ. More interestingly, we also obtain a quasi-isothermal tem- 
perature profile, which is a direct consequence of multi-phase cooling within such 
a potential well. 

1 Introduction 

Active cooling strongly regulates the evolution of the dynamically relaxed cep.- 
tral regions of many galaxy clusters, as inferred from X-ray observations 13. 
The analysis of X-ray spectra of these regions often exhibit an excess of low- 
energy photons, relative to a single-temperature spectrum. Recently, high 
EUV brightness excess have heen detected by EUVE pointed observations of 
the Virgo and Coma clusters a. Such observations strongly suggest the pres- 
ence of cold gas at a temperature well below the average X-ray temperature of 
the intra cluster gas. This may be an evidence for the presence of cold parcels 
of gas immersed in the hot, pervasive X-ray emitting gas. 

Gas cooling is a complex, unstable phenomenon. As a matter of fact, the 
gas is subject to the so-called cooling instability by which small temperature 
differences are amplified, leading to a clumpy structure on small scales. A long 
lived, tenuous hot phase can be maintained as the bulk of the mass cools down 
and eventually condense in cold clouds and, possibly, stars. This description 
of the thermal history of the gas has been investigated by amnerous authors, 
and is known has the "multi-phase cooling flow" formalismH'EKHD. 



In this paper, we intend to investigate multi-phase cooling using a fully 
hydrodynamical approach. We therefore proceed a step further than the usual 
stationary cooling flow approach_We describe the dynamical evolution of both 
gas and dark matter componentsH, starting from cosmologically relevant initial 
conditions, together with the thermodynamical evolution of the multi-phase 
medium, which ultimately leads to the condensation of cold clouds. We shall 
restrict ourselves to very high resolution, spherically symmetric simulations, 
since we are interested in the most central regions of X-ray clusters, where 
strong cooling takes place. 



2 Physical and Numerical Methods 

Nulsen (1986) derived the equations describing the evolution of the density 
distribution in a multi-phase medium. These multi-phase equations describe 
the flow at intermediate scales, in between microscopic scales, where atomic 
processes contribute to gas cooling, and macroscopic scales, relevant for the 
overall hydrodynamical evolution. Integrating the multi-phase equations over 
the density distribution, one obtains useful modified hydrodynamical equations. 
First, the continuity equation writes 

\± + V-u = ~(3(x,t) (1) 
p at 

where j3 is the mass deposition rate and p is the mean density in the 
macroscopic fluid element. Second, the energy equation writes 

1 dP 5„ 2 W ,n 2 A(T) 

PA^'^B^P (2) 

where P is the gas pressure, A(T) the cooling function evaluated at the mean 
temperature, and X(x,i) the cooling enhancement factor, due to the pres- 
ence of an underlying density distribution. Finally, the Euler equation writes 

= -VP - pV* (3) 

The usual single phase hydrodynamical equations are recovered with A = 1 
and (3 = 0. However, in the general multi-phase case, the functions 0(x, t) and 
X(x,t) have to be evaluated through a detailed treatment of the multi-phase 
distribution. This has been done so far only in the case of stationary cooling 
flowsLra. It has been however demonstrated that, given any reasonable initial 
multi-phase medium, a universal high-density cooling tail rapidly develops in 
the mass spectrum, leading to self-similar solutions □. In this case, the cooling 




Figure 1: Simulated profiles of a Coma-like cluster: enclosed overdensity, mean electron 
density, accretion rate and mean gas temperature at different redshifts (z=2.4, 1.2, 0.5, 0). 
Note that the X-ray temperature is here only 70% of the mean temperature. 



enhancement factor A is constant all over the flow, and the mass deposition 
rate f3 is exactly proportional to the mean cooling rate in the macroscopic fluid 
element. In this paper, we neglect the short relaxation phase anterior to this 
asymptotic cooling flow regime. We consider only the self-similar solutions 
found by Nulsen (1986), parameterized by a single parameter v. Details on 
our numerical treatment will be found elsewhere El. 



3 Results and Discussion 

We simulate the formation of a large X-ray cluster, similar to the Coma cluster 
(Mtot ^ 1 x 10 15 , f ~ 10 s ). We plot in figure (Q) the profiles we obtained at 
different redshifts. The dark matter density profile is similar to the Navarro, 
Frenk et White (1995) universal profile, with p oc r _1 in the center. We found 
that only large values of the multi-phase parameter are allowed (y > 3), in 
order to avoid a cooling catastrophe. This roughly corresponds to the limiting 
case v = +oo, which allows indeed a strict conservation of the gas mean 



entropy, as fluid elements sink towards the center. It turns out to be the 
only stable case. This case corresponds to a wide temperature distribution 
within the intracluster medium, with small parcels of gas of arbitrary high 
temperature. 

The X-ray emitting gas density profile show a typical core-halo structure, 
the core radius corresponding to the cooling radius of the cluster. Note that 
the central density decreases in time, due to continuous mass deposition in 
this strongly cooling region. This deposition process regulates the mass infall, 
and the central density remains at a value corresponding to t coo i ~ tn- The 
accretion rate steadily increases with radius as M oc r 2 , and reaches 2000 
Mq/ut in the outer regions, which is observed in strong cooling flow clusters. 
More interestingly, the mean temperature profile is nearly isothermal. This 
property is due, first, to the total mass profile (p oc which determines the 
gravitational potential, and, second, to the multi-phase model we considered 
(y = +00). We found also that the inner temperature profile is highly sensitive 
to the mass profile. This lead us to conclude that, within our multi-phase 
model, the total density profile p oc r _1 is the only one leading to an isothermal 
gas distribution. 

The multi-phase model (y = +00) we use here might be considered as 
unrealistic, since it involves a very wide temperature distribution in the intra- 
cluster gas. We believe on the contrary that this model mimics the continuous 
reinjection of high entropy gas by supernovae driven winds, a very important 
aspect of clusters thermodynamics. A more realistic treatment of multi-phase 
cooling is however necessary, and we intend to include within our hydrody- 
namical code a complete (multi-group) treatment of the mass deposition. 
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